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We use high order linked cluster series to investigate the hard core boson model on the triangular 
lattice, at zero temperature. Our expansions, in powers of the hopping parameter t, probe the 
spatially ordered 'solid' phase and the transition to a uniform superfluid phase. At the commensurate 
fillings n = |, | we locate a quantum phase transition point at (t/V) c — 0.208(1), in good agreement 
with recent Monte Carlo studies. At half-filling (n = |) we find evidence for a solid phase, which 
persists to t/V ^ 0.06. 

^ ' I. INTRODUCTION 
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There is considerable interest at present in the physics of a system of hard-core bosons on the two-dimensional 
triangular lattice, described by the Hamiltonian 



H = -ty^qjqj + atai) + V'^2'n l n j - pk^nj (1) 
O 1 to) to) 

o ; 

where each lattice site can be occupied by or 1 bosons (rij = 0,1). The first term represents hopping between 
nearest neighbour sites, the V-teiui is a nearest neighbour repulsion, n, = ajcij, and fj, is the chemical potential. This 
Hamiltonian also arises as a limiting case (U = oo) of the so-called 'Bose-Hubbard' model, where multiple occupancy 



^vO is allowed, albeit at a cost in energy. 

c — . 

Interest is this model has arisen, in the main, because of the possibility of an exotic 'supersolid' phase where 
long-range crystalline order and superfluidity may co-exist. Such a phase was conjectured over thirty year ago and 
has remained controversial and unobserved until, perhaps, very recently 1 . On the other hand evidence for supersolid 



phases has been found in lattice models of interacting hard-core bosons, and the possibility of experimental realizations 
with atoms in optical potentials 2 provides an exciting prospect of direct verification of the theoretical predictions. 
As is well known, a lattice model of hard-core bosons with nearest-neighbour repulsion can be mapped directly 
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to a spin-I antiferromagnet, in general in the presence of a magnetic field. Techniques developed and used for the 

Q , magnetic system can then be used to study the boson problem. The early work of Liu and Fisher 3 studied the phase 

^ . diagram of the model, on a bipartite lattice, using the usual mean-field approximation, and showed the existence of 
• i—i , 

a supersolid phase under certain conditions. More recently Murty et al. have used first-order spin-wave theory to 
study the model on the two-dimensional triangular and kagome lattices, and also find a stable supersolid phase. These 
lattices are attractive candidates for exotic quantum phases because of frustration and low dimensionality, and the 
consequent enhancement of quantum fluctuations. Both of these studies 3 ' 4 involve, however, approximations whose 
validity is difficult to judge. 

In an attempt to avoid uncontrolled approximations Boninsegni 5 used a Green Function Monte Carlo method to 
study the triangular lattice. He verified the presence of a supersolid phase, albeit within a smaller range of parameters 
than predicted at the mean-field level. In particular, he found for V/t > 7, a narrow supersolid phase at and above 
the commensurate filling n = h (n = (rii) is the average boson number density), and by particle-hole symmetry, also 
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at and below n = |. He found no supersolid phase at n = ^ for any value of V. More recent work 6,7,8,9 has further 
clarified the phase diagram, particularly at and near half-filling (n = ^). It appears, contrary to earlier expectations, 
that the supersolid phase persists at n = 5 although its spatial symmetry remains controversial. 

Our goal in this paper is to investigate the ground state properties of the triangular lattice hard-core boson model 
via the technique of high-order series expansions 10 ' 11 , an approach complementary to the Monte Carlo work. A 
particular advantage of the series approach is the ability to locate critical points reliably and accurately. To the best 
of our knowledge the hard-core boson problem has not previously been studied using this approach, although there 
has been some work on the Bosc-Hubbard model without nearest neighbour repulsion 12 . 

The paper is organized as follows. Section II discusses some general aspects of the model and briefly describes the 
series expansion methodology. Section III gives our results for the case of n = | and | , and discusses these in relation 
to existing work. Section IV considers the half- filling case n=\. Finally in Section V we summarize our results and 
attempt to draw conclusions. 

II. GENERAL FEATURES AND METHODOLOGY 

As already mentioned above, the hard-core boson Hamiltonian has a particle-hole symmetry. Under the transfor- 
mation di — of, fii = 1 — rii the Hamiltonian (1) becomes 

H = N(-zV - n) - t^2(dldj + d]di) + V ^ fiifij - (zV - n) ^ n t (2) 

(ij) (ij) i 

which, apart from a constant, is identical in form to the original. Here z(= 6) is the coordination number of the lattice. 
Thus the phase diagram is symmetric about the half- filled case n—\. Consequently our series for the commensurate 
filling n = | are unchanged for n = |, and we consider the two cases together in Section III. 
It is also useful to display the equivalent spin Hamiltonian, obtained via the transformation 

Oi^Sr. al^Sf, n l ^S z + \ (3) 

which yields, apart from a constant 

# sp i„ = V SiSj - ^(SfS^ + SfS]) - Ox - \zV) SI (4) 

(ij) (ij) » 

i.e., a Heisenberg model with exchange anisotropy (XXZ model) in a magnetic field. A phase with (S z ) ^ 0, (S x ) = 
then corresponds to a solid phase or normal fluid, while a phase with non-zero magnetization in the x-y plane 
corresponds to a superfluid (if (S z ) is uniform) or a supersolid (if (S z ) is spatially modulated). 
The series expansion method is based on writing the Hamiltonian in the form 

H = H + XH 1 (5) 

where H is exactly solvable and Hi is treated perturbatively. A linked-cluster approach 10,11 is the most efficient, 
permitting the series to be derived to quite high order. Here we take the nearest neighbour repulsion term as H , and 
the hopping term as the perturbation. The Hamiltonian (1) is particle conserving and hence the chemical potential 
term can be dropped in the calculations. For convenience we also set V — 1 to fix the energy scale. The perturbation 
parameter is then A = —t, and our series are obtained in powers of t. As the series are of finite length we are 
unable to effectively probe the large t region of the phase diagram. In magnetic language, our unperturbed state is 
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a commensurate Neel type state, and x-y spin fluctuations are included perturbatively. We are unable to treat the 
superfluid/supersolid phases directly. Nevertheless, as we shall show, a great deal of useful information about the 
model can be obtained. 
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FIG. 1: (a) The 3-sublattice decomposition of the triangular lattice; (b) Unperturbed configuration for n = | (the filled circle 
represent occupied sites); (c) Unperturbed configuration for n = |. 



The triangular lattice can be decomposed into three equivalent sublattices, which we denote a, (3, 7, as shown in 
Fig. 1(a). For filling n = | the unperturbed ground state has two sublattices (say a and (3) fully occupied and the 
third (7) empty This corresponds to a solid phase, with a density modulation with wavevector k = (4^,0), and is 
shown in Fig. 1(b). It is also referred to as the \/3 x \/3 phase. Although the ground state is threefold degenerate, 
these do not mix in any finite order of perturbation theory and we can start from one of the partners. Turning on 
the hopping term allows particles to occupy the other sublattice, and at some critical value t c the average sublattice 
occupancies will become equal. This represents the transition from solid to superfiuid. An order parameter can be 
defined as 

m = (n a ) — (n 7 ) = 3(n a ) — 2 (6) 

which equals 1(0) in the fully ordered (disordered) phase. The transition can also be located from the behaviour of 
nearest neighbour correlators C a p = {n a np), or from the static structure factor 

S(V = ^(n i n j )e i ^-^ (7) 

i,3 

at k = k*, where k* = (^,0) is the wavevector for the V3 x \/3 phase. For n = | a completely analogous treatment 
can be used, with the unperturbed ground state having one sublattice (say 7) fully occupied, and the other sublattices 
empty. 

At half- filling (n = |) the unperturbed system is equivalent to an isotropic Ising antiferromagnet in zero field, which 
is known to be disordered at T = 0, with a macroscopic ground state degeneracy. It seems, however, that a nonzero 
hopping term lifts this degeneracy and leads to a stable supersolid phase - an example of the 'order from disorder' 
phenomenon 13 , induced by quantum fluctuations. To study this case via series expansions we include an additional 
field in the unperturbed Hamiltonian (and subtract it via the perturbation term), to choose a unique ground state. 
We choose the striped state shown in Fig. 1(c), for which a two-sublattice decomposition of the triangular lattice is 
made. This will be discussed further in Section IV. 
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III. THE n= §,§ PHASES 

We turn now to our results, considering first the commensurate phase with n = |. The unperturbed ground state 
(•\/3 x V3) is shown in Figure 1(b). The order parameter m (Eqn. 6) will decrease with increasing t, and, assuming 
a normal 2nd order transition, will go to zero with a power law 

m ~ (t c - tf, t -» t~ (8) 

where (5 is a critical exponent. 

We have computed a series expansion for m, up to order i 12 . The series coefficients are given in Table 1. Dlog Padc 
approximants (shown in Table II) and integrated differential approximants 14 have been used to analyse this series, 
with the resulting estimates 

t c = 0.21(1), (3 = 0.09(1) (9) 

The estimate t c — 0.21 is consistent with the most recent Monte Carlo estimate 6 of 0.195 ± 0.025. We are unaware 
of any previous estimate of the exponent f3. The order and universality class of this transition remain unclear. Our 
exponent estimate ~ 0.09(1) is similar to that of the 3-state Potts model in two dimensions (/3 = 1/9). On the other 
hand the 3-state Potts model in three dimensions is believed to have a weak first-order transitions. 

It would be straightforward to compute the order parameter for specific values of t, and to display a curve of m(t), 
but we do not do this. 
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FIG. 2: The ground state energy per site Eq and (d 2 Eo/dt 2 ) 1 versus t for n = 2/3 filling, the results of several different 
integrated differential approximants to the series are shown. 

An alternative way of estimating the critical point is to compute the series for the ground state energy E a (t) and 
for the second derivative d 2 E /dt 2 (analogous to a 'specific heat'). This analysis is less precise, but the results are 
consistent with (9). In Figure 2 we show E and the inverse of the second derivative as functions of t, for < t < t c . 
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FIG. 3: The nearest neighbor correlators C a p and C Q7 and their difference, versus t, for n = 2/3 filling. Several different 
integrated differential approximants to the series are shown. 

The results of several different integrated differential approximants to the series are shown. The second derivative is 
seen to diverge at t <~ 0.21. 

We have also obtained series in t for the pair correlations Cij = (runj). One such series is given in Table I) - we are 
happy to provide others on request. Figure 3 shows three curves, as functions of t: the nearest neighbor correlators 
C a p and C Q7 , and their difference. The difference C a p — C a ~ ( decreases from 1 as t increases, and vanishes at the 
critical point where all sublattices are equally occupied. This gives an independent estimate of t c , consistent with (9). 

Another quantity of interest is the static structure factor 



This differs from the more usual definition (eq. Ref. 5) by the subtraction of the second term, which is necessary to 
implement the linked cluster method efficiently. The effect of this term is to remove the 5- function peak at k = 0, 
and reduce it at the special wavevector k*. At the transition point t Cl where (n^) is independent of sublatticc, the 
contribution of this second term vanishes, and our S(k) is expected to show a S- function peak at k*. This is confirmed 
in Figure 4, where we show <S'(k) curves along symmetry lines in the Brillouin zone for various t values. 

At k = k*, the Dlog Pade approximants to the series S^k*) give a critical point estimate t c = 0.208(3) with 
exponent —0.40(3) (see Table II). To obtain a more accurate estimate of 5(k*) versus t, we use integrated differential 
approximants to the series [5(k*)]~ 2 ' 5 (which vanishes linearly at t c ). The result is shown in Fig. 5, from which we 
can locate the critical point more accurately, at t c = 0.208(1). 
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FIG. 4: The static structure factor S(k) along high-symmetry cuts through the Brillouin zone for various t at n = 2/3 filling. 
The inset shows the Brillouin zone, and path chosen. 

IV. THE n= \ PHASE 

As noted in Section II, the case of half- filling (n = |) is rather more subtle, as the unperturbed system has an 
infinitely degenerate ground state. To allow a series expansion we introduce an auxiliary field to break this degeneracy. 
There are many ways to do this, but the simples choice is to decompose the triangular lattice into two sublattices 
(a, (3) shown in Figure 1(c), and to impose a staggered field, which favours occupancy of the a sublattice. Thus the 
unperturbed Hamiltonian is 

H = V^2ninj - h^^riirii (11) 

(ij) 1 

and the perturbation is 

Hi = -t^2(a\aj + a]oi) + h^^m (12) 

(ij) 1 

where rji = +1, —1 on sublattices a, (3 respectively. The full Hamiltonian is H = Hq + XV. Perturbation series are 
derived in power of A, and A = 1 corresponds to the original Hamiltonian. Note that the field terms cancel in this 
limit. The value of h is arbitrary and can be chosen to obtain optimal convergence of the series. 

We have derived series for the same quantities as in the previous Section: the ground state energy and the order 
parameter, m = (n a ) — (np), to order A 11 , pair correlations and the structure factor to order A 10 . The series are 
evaluated as single variable series in A, for fixed values of t, h. To estimate the critical value t c we proceed as follows. 
For any value of t the series will have a singularity at A = A c (t) and t c is determined from the condition A c (t c ) = 1. 
For t < t c , A c is larger than 1, implying that the system retains the striped order shown in Figure 1(c) at A = 1, 
while for t > t c , A c is less than 1. The analysis is less precise than for the n = |, | cases, but we obtain a consistent 
critical point at t c ~ 0.06. Table III presents a number of series, for the parameter choice t = 0.06, h = 0.25. Table 
IV presents the results of a Dlog Pade approximant analysis. The order parameter series, in particular, is quite well 
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FIG. 5: The static structure factor 5(k*) versus t at n — 2/3 filling. Several different integrated differential approximants to 
the series are shown. Since the series begins with a t 2 term we divide this out for convenience. 

behaved and shows a consistent sequence of estimates of A c converging to A c ~ 1.0. The series for correlators and for 
the structure factor give a similar conclusion, that t c ~ 0.06. 

In Figure 6 we show the structure factor along symmetry lines in the Brillouin zone for n = i. The ordering 
wavevector for the striped 2-sublattice ordered structure is k* = 27r(0, l/v 7 ^) and we see the development of a 
divergent peak at k*, as t — > 0.06. There is also a peak at the ordering wavevector for the V3 x V3 phase. Figure 7 
shows the divergence of S(k*) as t increases, towards 0.06. 

Our estimate for the transition t c <~ 0.06 is somewhat lower than the latest Monte Carlo estimate 6 of t c <~ 0.115. 
We discuss possible reasons for this in the final section of the paper. 

V. CONCLUSIONS 

We have used series expansions to investigate the ground state phase diagram of the hard-core boson model on the 
triangular lattice, which is of considerable current interest. Series are derived in powers of the hopping parameter t, 
and hence probe the spatially ordered phases and transitions to a uniform supcrfluid phase. Series have been computed 
for the ground state energy, the order parameters, the near neighbour pair correlations and the static structure factor. 
All of these quantities show consistent results with each other. 

At the commensurate densities n — |, | we locate a quantum phase transition at (t/V) c ~ 0.208(1), in good 
agreement with the most recent Monte Carlo results of 0.195 ± 0.025. The order parameter appears to vanish at 
(t/V) c with a power low, with an exponent (3 ~ 0.09. At n = \ (half-filling) early Monte Carlo work 5 indicated 
no ordered solid phase at any finite i, although the most recent work 6 ' 7,8,9 finds a supersolid phase persisting to 
t/V ~ 0.1. To derive series expansion for this case we stabilize one of the infinite manifold of ground states (a 
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'striped' 2-sublattice phase) by applying an external field. This phase, which presumably is supersolid, appears stable 
to t/V ~ 0.06. 

Using the present approach we are unable to probe off-diagonal long-range order and hence to study the superfluid 
phases. 

We are also unable to exclude the possibility of weakly first-order rather than second-order transitions. In this 
context we mention recent field theoretical work 15 , which argues that the transitions in this model lie in the class of 
'deconfined quantum critical points'. 
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FIG. 7: The static structure factor 5(k*) versus t at n = 1/2 filling. 
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TABLE I: Series for ground state energy per site Eo(t), order parameter m, difference of shortest correlator C a p — C a -y, and 
static structure factor S at k* for the n = 2/3 case. Nonzero coefficients t n up to order n = 12 are listed. 



n 


E (t) 


m 


Cap — Cot'f 


S(k') 





1.000000000 


1.000000000 


1.000000000 


0.000000000 


1 


0.000000000 


0.000000000 


0.000000000 


0.000000000 


2 


-1.000000000 


-2.250000000 


-2.750000000 


1.500000000 


3 


-1.000000000 


-4.500000000 


-5.500000000 


3.000000000 


4 


-9.500000000X10- 1 


-1.019750000X10 1 


-1.052083333X10 1 


8.168333333 


5 


-2.355555556 


-3.750666667X10 1 


-3.819305556X10 1 


3.518166667X10 1 


6 


-7.937379630 


-1.418886375xl0 2 


-1.498124375xl0 2 


1.428442162X10 2 


7 


-2.616389136X10 1 


-5.619014174xl0 2 


-5.933175977xl0 2 


6. 107654281 xlO 2 


8 


-9.005836612X10 1 


-2.307990196xl0 3 


-2.428296319xl0 3 


2. 667494106 xlO 3 


9 


-3.229654053 xlO 2 


-9.596323715xl0 3 


-1.008214369xl0 4 


1.170747249xl0 4 


10 


-1.201486606xl0 3 


-4.067064188xl0 4 


-4.264450526xl0 4 


5. 234712003 xlO 4 


11 


-4.607510926 xlO 3 


-1.750821492xl0 5 


-1.832564264xl0 5 




12 


-1.792870548xl0 4 


-7.572580451 xlO 5 
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TABLE II: [n/m] Dlog Pade approximants to the series of d 2 Eo/dt 2 , m, Cap — C al , and S(k*) for n = 2/3 case. The position 
of the singularity (pole), the critical index (index) from unbiased approximants are given. An asterisk denotes a defective 
approximant. 

n [(n-2)/n] [(n - l)/n] [n/n] [(n + l)/n] [(n + 2)/n] 

pole (index) pole (index) pole (index) pole (index) pole (index) 

d 2 E (t)/dt 2 

n=l 1.2500(-3.750) 0.0518(-0.006) 0.1834(-0.286) 0.1929(-0.350) 

n=2 0.2342(-0.388) 0.1514(-0.139) 0.1936(-0.357) 0.1953(-0.373) 0.1921(-0.345)* 

n=3 0.1726(-0.213) 0.1951(-0.371) 0.1933(-0.355)* 0.2124(-0.845) 0.2007(-0.438) 

n=4 0.2111(-0.691) 0.2020(^0.463) 0.2010(-0.444) 0.2036(-0.502) 

n=5 0.2006(-0.438) 0.2018(-0.459) 

m 

n=l 0.3333(0.5000) 0.2651(0.2517) 0.2138(0.1064) 

n=2 0.2750(0.2896) 0.1673(0.0261) 0.2209(0.1246) 0.2180(0.1162) 

n=3 0.2368(0.1741) 0.2257(0.1411) 0.2164(0.1109) 0.1999(0.0504) 0.2126(0.0986) 

n=4 0.1852(0.0262) 0.2103(0.0897) 0.2116(0.0946) 0.2108(0.0910) 0.2087(0.0809) 

n=5 0.2121(0.0965) 0.2111(0.0925) 0.2168(0.0994)* 0.2100(0.0876) 

n=6 0.2096(0.0860) 0.2104(0.0895) 

Cot/3 ~ Cotf 

n=l 0.3333(0.6111) 0.2884(0.3959) 0.2146(0.1213) 

n=2 0.2932(0.4219) 0.1049(0.0019) 0.2207(0.1393) 0.2182(0.1309) 

n=3 0.2431(0.2229) 0.2288(0.1712) 0.2166(0.1247) 0.1983(0.0487) 0.2134(0.1134) 

n=4 0.1906(0.0404) 0.2115(0.1047) 0.2123(0.1079) 0.2111(0.1025) 0.2090(0.0905) 
n=5 0.2124(0.1087) 0.2118(0.1058) 0.1988(0.0307)* 
n=6 0.2100(0.0970) 

son 

n=l 0.2902(-0.5805) 0.1508(-0.1568) 0.2228(-0.5051) 0.2001(-0.3291) 

n=2 0.1833(-0.2678) 0.1941(-0.3066) 0.2052(-0.3709) 0.2051(-0.3700) 0.2152(-0.5228) 

n=3 * 0.2051(-0.3700) 0.2052(-0.3709) 0.2072(-0.3895) 

n=4 0.2082(-0.4033) 0.2081(-0.4024) 
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TABLE III: Series for ground state energy per site Eo, order parameter m, difference of shortest correlator C aa — Cap and 
C aa — Cfjp, and static structure factor S at k* = (0, 27r/\/3) for t = 0.06 and the field h = 0.25 at n = 1/2 case. Nonzero 
coefficients A n up to order n = 11 are listed. 



n 


Eo 


m 


Caa C a fi 


Caa — Cfjf) 


5(k*) 





5.000000000 xl0~ 


i 


1.000000000 




1.000000000 


1.000000000 




0.000000000 




1 


0.000000000 




0.000000000 




0.000000000 


0.000000000 




0.000000000 




2 


-4.800000000x10" 


-3 


-1.280000000x10" 


-2 


-1.760000000xl0~ 2 


-1.280000000x10" 


-2 


1.280000000x10" 


-2 


3 


-1.984000000x10" 


-3 


-1.058133333x10" 


-2 


-1.454933333xl0~ 2 


-1.058133333x10" 


-2 


1.058133333x10" 


-2 


4 


-8.346453333x10" 


-4 


-6.818781867x10" 


-3 


-9.282071467xl0~ 3 


-6.818781867x10" 


-3 


6.969821867x10" 


-3 


5 


-3.763032178x10" 


-4 


-4.322561214x10" 


-3 


-5.768006277xl0~ 3 


-4.322561214x10" 


-3 


4.758822760x10" 


-3 


6 


-1.920603977x10" 


-4 


-2.973745768x10" 


-3 


-3.879098112xl0~ 3 


-2.973745768x10" 


-3 


3.673001969x10" 


-3 


7 


-1.129137214x10" 


-4 


-2.252261754x10" 


-3 


-2.888143560xl0~ 3 


-2.252261754x10" 


-3 


3.111869445x10" 


-3 


8 


-7.425851217x10" 


-5 


-1.819385962x10" 


-3 


-2.309557263xl0~ 3 


-1.819385962x10" 


-3 


2.744159087x10" 


-3 


9 


-5.242495609x10" 


-5 


-1.520557155x10" 


-3 


-1.918242621xl0~ 3 


-1.520557155x10" 


-3 


2.454367694x10" 


-3 


10 


-3.867629129x10" 


-5 


-1.296181354x10" 


-3 


-1.627204515xl0~ 3 


-1.296181354x10" 


-3 


2.215862063x10" 


-3 


11 


-2.946020021x10" 


-5 


-1.122501973x10" 


-3 
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TABLE IV: [n/m] Dlog Pade approximants to the series of d 2 Eo/d\ 2 , m, C aa — C a p, C aa — Cpp, and S(k*) for t = 0.06 and 
the field h — 0.25 at n — 1/2 case. The position of the singularity (pole), the critical index (index) from unbiased approximants 
are given. An asterisk denotes a defective approximant. 



11 


[(" - 2)/n] 


[(n-l)/n] 


[n/n] 


[(n + l)/n] 


[(n + 2)/n] 




pole (index) 


pole (index) 


pole (index) 


pole (index) 


pole (index) 








d 2 E /d\ 2 






n=l 




2.2586(-2.8007) 


1.4547(-1.1617) 


0.9731(-0.3478) 


0.9220(-0.2803) 


n=2 


1.6180(-1.5630) 


0.7556(-0.1148) 


0.9164(-0.2716) 


0.9458(-0.3145) 


1.1581(-6.1720)* 


n=3 


0.9303(-0.2941) 


1.0207(-0.4982) 


1.0060(-0.4512) 


0.9990(-0.4284) 


0.9929(-0.4076) 


n=4 


1.0063(-0.4522) 


0.9889(-0.3895) 


0.9707(-0.3098) 






n=5 


0.9711(-0.3113) 










n=l 




* 


m 

0.8065(0.0166) 


1.1500(0.0483) 


1.2383(0.0649) 


n=2 


* 


1.3497(-0.0279)* 


1.2735(0.0755) 


1.2075(0.0577) 


1.3029(0.0748)* 


n=3 


1.1914(0.0521) 


1.1288(0.0392) 


1.0916(0.0316) 


1.0711(0.0271) 


1.0576(0.0241) 


n=4 


1.0746(0.0277) 


1.0591(0.0243) 


1.0364(0.0187) 


1.0094(0.0124) 


0.9977(0.0100) 


n=5 


0.9887(0.0083) 


0.9929(0.0090) 


0.9946(0.0094) 






n=6 


0.9962(0.0097) 
















Caa Ca(3 






n=l 




* 


0.8065(0.0229) 


1.1563(0.0675) 


1.2532(0.0931) 


n=2 


* 


1.3327(-0.0378)* 


1.2969(0.1123) 


1.2229(0.0830) 


1.3011(0.1043)* 


n=3 


1.2075(0.0751) 


1.1419(0.0559) 


1.1015(0.0443) 


1.0783(0.0373) 


1.0631(0.0326) 


n=4 


1.0820(0.0382) 


1.0643(0.0327) 


1.0407(0.0251) 


1.0149(0.0171) 




n=5 


1.0014(0.0133) 


0.9989(0.0127) 














Caa — C/3/3 






n=l 




* 


0.8065(0.0166) 


1.1500(0.0483) 


1.2383(0.0649) 


n=2 


* 


1.3497(^0.0279)* 


1.2735(0.0755) 


1.2075(0.0577) 


1.3029(0.0749)* 


n=3 


1.1914(0.0521) 


1.1288(0.0392) 


1.0916(0.0316) 


1.0711(0.0271) 


1.0576(0.0241) 


n=4 


1.0746(0.0277) 


1.0591(0.0242) 


1.0364(0.0187) 


1.0094(0.0124) 




n=5 


0.9887(0.0083) 


0.9929(0.0090) 














5(k*) 






n=l 




2.0378(-1.6846) 


1.2297(-0.6135) 


0.9509(-0.2837) 


0.9605(-0.2952) 


n=2 


1.4031(-0.8844) 


0.8620(-0.1860) 


0.9602(-0.2948) 


0.9482(^0.2815)* 


1.0447(-0.4819) 


n=3 


0.9726(-0.3144) 


1.0375(-0.4476) 


1.0246(-0.4121) 


1.0163(-0.3890) 




n=4 


1.0250(-0.4132) 


0.9918(-0.3047) 









